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1. INTRODUCTION 

An island with these loads is formed by a decentralized generator that is separated from the main 
power grid but continues to supply the loads. This is referred to as islanding. This can be deliberate or 
unintentional. Indeed, the disconnecting of the electrical network during a repair operation can result in the 
islanding of generators. The island is recognized and can be switched off by intervening individuals because 
the network loss is voluntary. 

The inadvertent islanding, which was caused by a power outage, is of higher interest. It does, in fact, 
pose a number of risks, including the possibility of electrical equipment damage due to large voltage and 
frequency variations, difficulties in phase shifting between the grid and the distributed generator when 
reconnecting after a fault, and risks to maintenance personnel and the general public. It is consequently 
critical to detect any islanding situations and lower the islanded system's operational time. 

Three techniques have been proposed in this article. Their operation is based on the following 
characteristics: frequency, amplitude or phase of the voltage at the grid connection point of the decentralized 
generator. This study proposes an inventory of the different methods of islanding, and highlights their 
detection principles, advantages and disadvantages. To test the performance and reliability of these 
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techniques, as well as to check the performance and efficiency of the effective method plan, simulation in the 
physical security information management (PSIM) environment was established. 


2. ACTIVEISLANDING DETECTION METHOD 
2.1. Active frequency drift (AFD) 

With active frequency drift (AFD), the inverter current waveform is very slightly unstable, as shown 
in Figure 1, for a drift regime. AFD can be obtained by including a zero t; < 833s (dead time of PV 
inverter current) bending segment. Or by forcing the frequency of the current to be higher than the frequency 
of the voltage of the previous cycle by a frequency offset df < 1.5Hz. The value of the zero current segment 
t, or the frequency dfis limited to the total harmonic distortion of the current THD; < 5%. When there is a 
network, it will maintain the voltage frequency point of common coupling (PCC). However, when the 
photovoltaic system is disconnected from the grid [1], the voltage frequency at PCC tends to drift upwards, 
reaching a value higher than the resonance frequency fo the parallel resistor-inductor-capacitor (RLC) load 
[2]. This increases the risk that the net inductive load will trip the over freuency protection (OFP) at the 
nominal grid frequency. In both cases, the basic component of the current flows through a small angle Oarp 
(Equivalent phase angle of current and voltage of inverter with AFD) before the voltage. 

The waveform of the distorted current carries out the active frequency drift upwards, as shown in 
Figure 1. Displays the PCC voltage, the waveform of the inverter output current, and the basic components of 
the inverter current. The dominant equation is (1): 


ipy = Isin(2rf;t) = Isin[2n(f + df)]t (1) 


where f = f, is the frequency of the PCC voltage. 
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Figure 1. Current and voltage waveforms of upward AFD [2] 


Figure 1 shows that although the output current of the photovoltaic (PV) inverter is in phase with the 
PCC voltage, the waveform of the current is neither odd nor even. In order to simplify the calculation of the 
phase angle between the basic components of voltage and current, the beginning of the basic current has been 
shifted by—t,/2. Therefore, the basic component of the output current of the ipy inverter drives the PCC 
voltage at an angle proportional to —t,/2 [3], [4]. The angle can be calculated in (2): 


tz 


2. _ AFD 
Ty 2m (2) 
where t; = 2a es *_and is commonly called dead time T,is the time period of the voltage. 
fo fi f fréf 


Therefore, the phase angle between the fundamental current of the inverter and the PCC voltage is (3). 


5; _ 11 _ not 
barp = nft, = nf G aa) = fief (3) 


This phase angle plays an important role in mapping the non-detection zone (NDZ). AFD can 


achieve acceleration or penetration properties in theory. Their impact, however, is highly dependent on the 
sort of reactive stress. The load voltage, or PCC voltage, is specified by the grid voltage before disconnecting 
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from the grid, and the inverter current is constrained to be in phase with the PCC voltage. The current will be 
distorted due to the predetermined acceleration. The load voltage will lose the reference voltage after being 
disconnected from the grid and will no longer be fixed by the grid. If the load is a pure resistance, the preset 
acceleration rate will raise the frequency of the inverter current in the next cycle, while also raising the 
frequency of the load voltage. The voltage waveform will lag behind the inverter current because the load is 
capacitive, and the effective acceleration will be smaller than the specified value. As a result, the generated 
acceleration increment is smaller than the preset value, and as the frequency increases, it gets smaller and 
smaller. They can reach 0 before the trip limit is reached, resulting in islanding. The effective frequency will 
increase for net inductive loads, which will assist speed up the process and improve islanding protection. 
Because most loads are inductive, AFD's preferred mode of operation is acceleration [5], [6]. 


2.2. Slip mode phase shift (SMS) 

The difference between the frequency of the PCC voltage in the previous cycle and the grid's 
nominal operating frequency drives the change in current phase angle in the slip mode phase shift (SMS) 
approach. When the frequency increases, phase advance happens, and when the frequency decreases, phase 
delay occurs [7], [8]. 

The refernce of the inverter current is (4): 


ipv = Isin(2rft + sms) (4) 


with, 9545: With SMS, the current and voltage phase angles of the inverter are equivalent. Set the phase shift 
as a sinusoidal function of the frequency deviation from the mains frequency fg, then: 


n f-fg ) (5) 


27 . 
Osus = zeg Omsin (z = 


with, Om: Maximum phase shift SMS (°). Where fmis the frequency at which the maximum phase shift 
fmoccurs, and f,,is the maximum phase shift in degrees (°). Like AFD, this phase angle plays an important 
role in NDZ mapping. 

After connecting the utility program, it will provide the solid phase and frequency, which are 
considered as the constant voltage waveform reference of the PCC. However, for an island system, the 
frequency at the PCC will vary according to the current injected by the inverter and will reach a steady state 
value, which can be calculated according to the load and SMS characteristics of the inverter. The steady-state 
frequency can be calculated as the frequency at which the inverter angle and the load have the same value. 
Frequency can also be identified by graphs [9], [10]. the relationship curve between the load angle roaa and 
frequency intersects the curve between the SMS angle 95,5 and frequency according to the formula. Figure 2 
shows the frequency response SMS and the load phase response as a function of frequency. This figure 
shows the phase response of two loads with the same resonant frequency, where fọ = 60Hz and the quality 
factors are Q; = 5 and Q; = 1.5. The intersection of the SMS group and the load phase response is indicated 
by points A, B, and C [11], [12]. 
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Figure 2. SMS phase response cuvre and parallel RLC load phase response cuvre [12] 


However, the intersection of the load line and the phase response box of the inverter may be a stable 


i ; . do dé f í 
or unstable operating point. The unstable operating point of re > Ea is at the point A where the load 


quality factor is Q; = 1.5, and the stable operating point of is at the point A where the load quality factor is 
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Qr = 5. Points B and C are stable with respect to the load Q; = 1.5. The stable operating point indicates the 
steady state frequency of the island system. The SMS parameters must be selected in such a way that there is 
no stable operating point within the permissible range of the mains frequency between the settings of the 
under frequency protection (UFP) and OFP devices. 


2.3. Sandia frequency shift (SFS) 

The current approach of islanding detection is a continuation of the frequency bias method, also 
known as active frequency drift with positive feedback. Positive feedback is used in the Sandia frequency 
shift (SFS) approach. This feedback will be applied to the voltage frequency of the interconnection point, and 
will try to generate interference on the grid by changing the frequency of the inverter output current to detect 
islanding more quickly. In order to set the inverter frequency to a different strategy than the grid frequency, 
as shown in Figure 3, the inverter output current shape is smoothly added or switched off [11], [12]. 

Figure 3, shows the example of the PV system inverter with an undistorted sine waveform for 
comparison. During the first half of the cycle, the output current of the PV system is a sinusoidal waveform 
whose frequency is higher than the frequency of the PCC voltage. When the output current of the PV inverter 
reaches zero, it remains at zero for the time t, before the start of the second half of the cycle. The higher 
frequency also applies to the second half of the current cycle. The dead time of the second cycles is not set by 
the PV inverter [13], [14]. 
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Figure 3. Waveforms with the SFS method [11] 


A positive reaction is used to increase cutting factor cy, which is mathematically defined as the ratio 


between time zero t, and half of the period of the PCC voltage waveform =. Therefore, 
2tz 
cp = (6) 


with an increasing deviation of the frequency from the nominal value. The erroris generally selected to be a 
linear function of the frequency of the increasingthe PCC voltage. 


Cp = Cy + F(Af) = Co + keps(f — fa) (7) 


Where Cro is a chopping factor when there is no frequency error, and k,r, is an accelerated gain 
[15], [16]. The SFS waveform has neither odd nor even symmetry, consequently, there will be a phase shift 
of the fundamental component of this form, wave equal to 0.5wt. If the zero current segments are small, the 
componentsharmonics of the current are also small, the SFS current can beapproximated by its set of out-of- 
phase fundamental components. Thus, the angle of the fundamental component of the inverter current varies 
with the frequency of the PCC, and the cutting factor cf [17], [18]. 


mcr(f) 
2 


wtz 


Osrs(f) = -~ = aft = 


(8) 


In order to remain in phase with the network voltage, at each zero crossing of the voltage of the grid 
the waveform of the inverter current is reset [19], [20]. Therefore, when the frequency of the current is higher 
than the mains frequency, the waveformof the current will complete its half cycle before the mains voltage 
drops to zero. Therefore, the current waveform will remain at zero until the next zero crossing ofthe mains 
voltage, then it will start a new sine half cycle [21]. Byelsewhere, if the frequency of the current is lower than 
that of the mains voltage, the passage at zero voltage will precede the end of the current half-cycle. Thus, this 
will resetthe current and will start a new sinusoidal cycle [22], [23]. 
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Therefore, this will produce acut-off in the waveform of the inverter output current. As shown in 
Figure 4. Currently, this method appears to be one of the most proactive methods ofeffective in preventing 
islanding. In fact, it has one of the weakest NDAs. Frommore, when several inverters are connected in 
parallel to the same node, it is easyto be installed, inexpensive and always effective [24]. On the other hand, 
if the system is connected toa weak grid where the source of distributed generation is of some importance, 
thismethod will result in a slight decrease in the quality of the waveform of the current fromoutput of the 
inverter and can cause undesirable energy instability. The network's main one [25]. Therefore, this method 
remains a good compromise between the quality of the waveform, the efficiency of islanding detection and the 
effect on the transient response of the system. However when the local load connected to the point of 
interconnection has a factor ofhigher quality, the effectiveness of this method will be considerably reduced [26]. 


3. SIMULATION RESULTS AND DISCUSSION 

A simulation was carried out for different detection methods with different loads in order to verify 
the theoretical analysis for the NDZ. It should be noted that if islanding is to occur, the voltage of the PCC 
must remain within the voltage limits when the network is disconnected. Thus, the amplitude of the reference 
current of the PV inverter is defined by ipymax = V2V, /R for the simulation. AFD performance is 
implemented for the case of frequency drift 6f = 1 Hz, and SMS is implemented for the case ofêm = 10°, 
fm = fg = 3 Hz and SFS is implemented forcfo = 0.05, ksfs = 0.05. 


3.1. Simulation result of the AFD method 

Figure 4, shows the time response of the system with the AFD method, the PCC voltage, the current 
inp the fault signal and the frequency fpcc, when the circuit breaker of the network is open and the local load 
is R = 14.4, C = 460.52 uF, and L = 15.28 mH, which gives the resonance frequency fọ = 60 Hz and the 
quality factorQ, = 2.5. The mains circuit breaker is open, i.e. the mains is disconnected at t = 0.07083 s. 
From the voltage waveform it can be seen that with the preset reference value of the current ipymax = 
V2, / R = 11.79 A, the voltage amplitude has not had any apparent change before and after the mains 
disconnection. However, the frequency of the voltage fpcc has increased and exceeded the upper limit of the 
second cycle frequency after disconnection from the grid. It should be mentioned that the indicates the 
frequency of the last voltage cycle. The PV inverter stopped feeding the local load within six cycles after the 
frequency exceeded the upper limits of t = 0.1822 A by switching off the inverter switches. It can be seen that 
the current ipy and the voltage Vpcc all dropped to zero after the fault signal was changed to 1. Thus, a load 
with f0 = 60 Hz and Q; = 2.5 is outside of the AFD NDZ during ôf = 1 Hz, which makes it possible to 
check the theoretical result. 
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Figure 4. Simulation result of AFD supplying a local load Pioaa = 1 kW, with fo = 60 Hz and Qp = 2.5 


In order to verify that the ZDN covers all the different resistance values, another simulation is 
performed with the load R=28.8, C = 187.33 uF, and L=38.85 mH, which gives the resonance frequency 
fo = 59 Hz and the quality factor Q; = 2. Figure 5 shows the time response of the system, the PCC voltage, 
the current ipp, the fault signal and the frequency fpcc. The mains circuit breaker is open, i.e. the mains is 
disconnected, at t=0.07083 s. The amplitude of the reference current is set to Ipym, = 5.89 A to maintain the 
amplitude of the PCC voltage. The frequency of the voltage fpcc has decreased to a steady state value within 
the frequency protection threshold. The fault signal remains at 0. The undetectable islanding operation has 


occurred, and the inverter is continuously feeding the load, resulting in the PV current i pv and the Vpcc 
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voltage. Thus, a load with fọ = 59 Hz and Qp = 2 is in the AFD NDZ for f = 1 Hz as predicted by the 
theoretical analysis. 


Figure 5. Simulation result of AFD supplying a local load Pioaq = 0.5 kW, with fo = 60 Hz and Qp = 2 


3.2. Simulation result of the SMS method 

Similar tests were performed for the grid-connected inverter using SMS with Om = 10° and 
fm = fg = 3 Hz. Figure 6, shows the system response in the time domain with SMS, PCC voltage, current 
ipy» fault signal and frequency fpcc, when the grid breaker opens and the local load is R = 14.4, 
C=460.52 uF and L=15.28 mH, which gives the resonance frequency fọ = 60 Hz and the quality factor 
Q; = 2.5. The mains circuit breaker is open, i.e. the mains is disconnected, at t = 0.07083 s, the positive 
peak of the fifth cycle of the PCC voltage. The frequency of the voltage fpcc has decreased and has become 
lower than the setpoint UFP after disconnection from the grid. The PV inverter stopped feeding the local load 
within 6 cycles after the frequency exceeded the lower frequency limits at t=0.4027 s by turning off the 
inverter switches. It can be seen that the PV current and the Vpccvoltage all dropped to zero after the fault 
signal changed to 1. Thus, a load with fọ = 60 Hz and Q; = 2.5 is outside the NDZ SMS for 6,, = 10°, so 
the results of the calculation can be checked. 
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Figure 6. Simulation result of SMS supplying a local load Pioaa = 1 KW, withfo = 60 Hz and Qp = 2.5 
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Another simulation is performed with the load, C = 552.62 uF and L = 12.73 mH, which gives the 
resonance frequency fo = 60 Hz and the quality factor Q; = 3. Figure 7 shows the time domain response of 


the PCC system voltage, ip, current, fault signal and frequency fpcc. The mains circuit breaker is open, i.e. 


the mains is disconnected at t = 0.07083. After disconnection from the grid, the frequency of the fpcc voltage 
has decreased, but has reached the steady state within the frequency protection threshold. The fault signal 
remains at 0. The undetectable islanding operation has occurred, and the inverter is continuously feeding the 


load, as shown by the current ip, and the VPCC voltage. Thus, a load with fọ = 60 Hz and Qr = 3 and is 
within the NZD range of the SMS for m = 10° as predicted by the theoretical analysis. 
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Figure 7. Simulation result of SMS supplying a local load Pioaa = 1 kW, with fo = 60 Hz and Q; = 3 


3.3. Simulation result of the SFS method 
Similar tests were performed for the grid-connected inverter using FSS with cfo = 0.05, 


ksrs = 0.05, Figure 8, shows the time domain response of the system with FSS, PCC voltage, ipy current, 
fault signal and frequency fpcc. When the grid breaker opens and the local load is R = 14.4, C = 460.52 uF 
and L = 15.28 mH, which gives the resonance frequency fọ = 60 Hz and the quality factor Q; = 2.5. The 
circuit breaker of the open network, i.e. the network is disconnected, at t = 0.07083 s, the positive peak of the 
fifth cycles of the Vpcc voltage. The frequency fpcc of the voltage has increased and exceeded the lower 


limit of the frequency after disconnection from the grid. The PV inverter stopped feeding the local load 
within 6 cycles after disconnection from the grid. The frequency exceeded the upper frequency limit at 


t = 0.1626 s by switching off the inverter switches. It can be seen that the PV current ipy and the Vpcc 


voltage all decreased from zero after the fault signal changed to 1. Thus, a load with fo = 60 Hz and 
Q; = 2.5 is outside the FSS NDZ area for Cfo = 0.05 and ksrs = 0.05, which allows checking the results of 
the calculation. 
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Figure 8. Simulation result of SFS supplying a local load Poaa = 1 kW, with fọ = 60 Hz, and Qp = 2.5 


Another simulation is performed with the load, C = 561.04 uF and L = 12.93 mH, which gives the 
resonance frequency fọ = 59.1 Hz and the quality factor Q; = 3. The time domain response of the PCC 


system voltage, ipy current, fault signal and frequency fpcc is shown in Figure 9. The mains circuit breaker 


is open, i.e. the mains is disconnected at t = 0.07083 s. The frequency of the voltage f PCC has reduced after 
being disconnected from the grid, but it had now stabilized within the frequency protection level. The fault 
signal is still at zero. The undetectable islanding operation has occurred, and the inverter is continuously 


feeding the load, as shown by the current ipy and the Vide voltage. Thus, a load with fọ = 59.1 Hz and 


Q; = 3is in the NDZ zone of the SMS for Cfo = 0.05 and ksrs = 0.05 as predicted by the theoretical 
analysis. 

One technique was picked after an investigation of the best islanding detection methods. The 
method chosen is actually a combination of a passive method and two active methods. The option of a 
combination of methods was kept since decreasing the NDZ improves the performance of the islanding 
detecting system. 
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Figure 9. Simulation result of SFS supplying a local load Pioaa = 1 kW, with fo = 59.1 Hz and Q; = 3 


The system's basic protection will be provided via the under/over voltage and under/over frequency 
methods. Because the frequency or amplitude of the voltage at the point of connecting exceeds the prescribed 
limitations, the inverter stops sending electricity to the grid. The root mean square (RMS) value of the 
voltage amplitude at the point of interconnection must be between 88% and 110% of the nominal grid 
voltage, and the frequency must be between 59.5 Hz and 60.5 Hz, according to the European standard. 

The active SFS approach was chosen as the second method. This strategy was chosen because it 
appears to be one of the most effective active strategies for avoiding islanding scenarios currently available. 
It does, after all, have one of the lowest NDZs. When many inverters are paralleled on a single connection 
node, it is also simple to construct, economical, and effective. 


4. CONCLUSION 

According to their placement in distributed generation systems, these approaches can be split into 
two categories: local and remote procedures. The distant sensing technique is used on the main grid side, 
whereas the sensing algorithm is used in the inverter for local techniques. Furthermore, local approaches are 
divided into two categories: i) passive methods, which rely on parameter estimation such as voltage, current, 
and harmonics, and ii) active methods, which use disturbances to detect islanding. Simulation was used to 
validate the theoretical analysis and the equation characterizing the NDZ of three common active integrated 
database management systems (IDMs). Matlab/Simulink was used to model and simulate the grid-connected 
PV system. The hysteresis control blocks, frequency and RMS, and the reference current with IDM are all 
model components. To verify the theoretical study, simulations of various IDMs, including AFD, SMS, and 
FSS, were run with various resonant frequencies and quality factor loads. These three IDMs' NDZ simulation 
results were plotted to compare with the theoretical results. The comparison of the various simulation results 
demonstrates the efficiency of the SFS method over the other methods. Because it has one of the lowest 
NDZs, this technique was chosen. 
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